Momentum distribution and correlation of two-nucleon relative motion in 6 He and 6 Li 
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The momentum distribution of relative motion between two nucleons gives information on the 
correlation in nuclei. The momentum distribution is calculated for both 6 He and 6 Li which are 
described in a three-body model of a+N+N. The ground state solution for the three-body Hamil- 
tonian is obtained accurately using correlated basis functions. The momentum distribution depends 
on the potential model for the N-N interaction. With use of a realistic potential, the 6 He momen- 
tum distribution exhibits a dip around 2fm~ 1 characteristic of 5- wave motion. In contrast to this, 
the 6 Li momentum distribution is very similar to that of the deuteron; no dip appears because it is 
filled with the D-wave component arising from the tensor force. 

PACS numbers: 27.20.+n, 21.60.-n, 25.60.Gc, 21.10.Gv 



I. INTRODUCTION 

In this article we study the momentum distribution 
and the correlation of relative motion of the halo neu- 
trons in two-neutron halo nuclei. The importance of 
the correlation in nuclei is widely recognized; obviously 
the correlation plays a vital role in binding a Borromean 
three-body system which has no pairwise bound states. 
Recently an experiment using the technique of intensity 
interferometry [l| has been done in order to extract the 
spatial correlation function of the two neutrons in the 
halo nuclei such as 6 He, n Li and 14 Be. In this exper- 
iment the momenta of the two neutrons and the core 
nucleus are measured after the dissociation of the halo 
nucleus. 

Another type of experiment has very recently been per- 
formed at RIKEN to probe the spatial correlation in 6 He 
and 6 Li, and the analysis of data is in progress Q. This 
experiment utilizes the one-nucleon exchange process in 
the interaction of these A=6 nuclei with a proton. The 
basic idea is to utilize the well-established one-nucleon 
exchange process which is observed at the backward angle 
in the proton-deuteron elastic scattering p|. Expecting 
this mechanism to occur in the A~6 nuclei, the cross sec- 
tion for the reaction 6 He(p, dn)a ( e Li(p,dp)a) has been 
measured under the backward scattering kinematics of 
the proton and two-nucleon system in 6 He ( 6 Li). Un- 
der the assumption of a quasi-elastic approximation, the 
cross section is expected to be sensitive to the relative 
momentum between the two nucleons Q. Though the 
reaction mechanism may not be as simple as expected, a 
theoretical analysis of the momentum distribution in 6 He 
and 6 Li should be important as a first step to understand 
the physics involved in the experiment. 

There are a number of calculations for the structure of 
6 He in various models. We employ a three-body model 
of a + N + N for 6 He and 6 Li. It will be important for 



our purpose to use a realistic potential for the interaction 
between the valence nucleons because the correlation be- 
tween them is primarily determined by the N-N potential 
and the distribution at high momentum should be sensi- 
tive to the short-ranged repulsion and tensor force. We 
refer to Refs. 0,0] as the most relevant calculations which 
include the tensor force. However, none of the literatures 
including them has considered the momentum distribu- 
tion between the valence nucleons. Of course there are 
many studies which have investigated the momentum dis- 
tribution of fragments such as the core nucleus or the 
nucleon [7|. 
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The plan of this paper is as follows. In Sec. [TT| we 
present our three-body model and some of the details 
on explicitly correlated basis functions which are used to 
solve the three-body problem. Some formulas needed to 
compare 6 He and 6 Li properties are given in this section 
with emphasis on the two-nucleon correlation function 
and the momentum distribution. The latter is best de- 
fined through the Wigner distribution function 8]. One 
well-known formulation of single-nucleon momentum dis- 
tribution near nuclear surface was performed by Hufher 
and Nemes @ using the Wigner function. Results of 
calculation are presented in Sec. IIHI together with some 
discussions. Those included are the spectroscopic prop- 
erties, the two-nucleon correlation functions which re- 
veal "dineutron" as well as "cigar-like" configurations 
and the comparison of the momentum distributions. We 
will show that the tensor force which works differently 
between 6 He and 6 Li plays a key role in the momen- 
tum distribution around 2fm . A conclusion is given 
in Sec. HVl In Appendix we give a formula to calculate 
a density matrix or the Wigner function using the corre- 
lated basis functions. 
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II. FORMULATION 
A. Three-body Hamiltonian 

The wave functions for 6 He and 6 Li are deter- 
mined from variational calculations for the core (a 
particle) +N+N three-body system which is specified by 
the Hamiltonian 

H = T r +T R + U 1 + U 2 + v 12 . (1) 

The subscripts of the kinetic energies T stand for the 
relative distance vector r between the two nucleons, and 
the relative distance vector R from the a particle to the 
center of mass of the two nucleons. The set of Jacobi 
coordinates (r, R) is called T-type hereafter. The poten- 
tial Ui is the N-a potential and V12 is the N-N potential. 
The a particle is treated as a structureless particle. 

As the two-nucleon potential Vi2, we use a realistic 
potential, G3RS (Gaussian soft core potential with 3 
ranges) case 1 potential [l0|], which contains central (C), 
spin-orbit (LS) and tensor (T) terms. In fact, the G3RS 
potential contains L 2 and quadratic L ■ S terms in even 
L waves. Their contribution is small; the deuteron en- 
ergy with (without) these terms is —2.17 (—2.28) MeV. 
Thus we ignore these L-dependent terms in what follows. 
The D state probability of the deuteron is 4.8% in the 
G3RS potential. We use the G3RS potential because 
its Gaussian radial form makes a numerical calculation 
much faster than, e.g., the AV8' potential [llj. It is in- 
structive for the study of the N-N correlated motion to 
compare results of calculation between the realistic po- 
tential model and an effective potential model which has 
a mild short-ranged repulsion. As such an effective po- 
tential we employ the Minnesota potential [l^| . abbre- 
viated to MN, which has no tensor component. This 
potential renormalizes the effect of the tensor force into 
the central force and reproduces the binding energy and 
the root mean square (rms) radius of the deuteron. 

As for the N-a potential Ui, we adopt a phenomeno- 
logical potential [13[, abbreviated to KKNN potential, 
which is determined so as to simulate the nonlocal po- 
tential that derives from a microscopic calculation in an 
N+a cluster model. The KKNN potential is parity- 
dependent, contains the central and spin-orbit compo- 
nents, and reproduces very well the low-energy N-a scat- 
tering phase shifts of S and P waves. The potential 
is slightly too repulsive in D and F waves [14j. The 
Coulomb potential for p-a is taken as 

^ Co >) = v erf (^ r )' w 

with =0.257 far 2 . The KKNN potential is deep 
enough to support an S'-wave bound state; the state is 
considered redundant and must be removed in the three- 
body calculation because no bound states exist for 5 He 
and 5 Li. 



B. Variational solution with correlated Gaussians 

Trial wave functions for the ground states of 6 He and 
6 Li are expressed, in LS coupling scheme, as a combina- 
tion of explicitly correlated Gaussians: 

K 

*jm =^C i #j M (A i , J 4 i ,K i ), (3) 
i=i 

with the basis function 
^jm(A={LS),A,u,x) 

= (1 - P12) {e-i* A * [y L (ux) X s(l, 2)U Vtm t (1, 2)} . 

(4) 

Here the permutation P12 ensures the antisymmetry of 
the valence nucleons. 

The basis function is specified by a set of nonlinear pa- 
rameters, the orbital and spin angular momenta A=(LS), 
a 2x2 positive-definite, symmetric matrix A, and a 2x1 
matrix u. The symbol ~ indicates the transpose of a ma- 
trix, and the square bracket [J-^Xs] denotes the angu- 
lar momentum coupling. The short-hand notation xAx 
with x=(xi,X2) stands for A\\x\ J r2A\2X\ ■ X2+A22X 2 ,, 
where the coordinates X\ and X2, called V-type, are the 
distance vectors of the valence nucleons from the a parti- 
cle; X\=R+\r and X2=R— \r. The exponential part of 
the basis function is rotation-invariant. The cross term 
A\2X\-X2 describes explicitly the two-nucleon correlation, 
which is vital to obtain a precise solution in a relatively 
small basis dimension [l5| . The angular part of the basis 
function is expressed by the solid spherical harmonics, 
yLM L {ux) — \ux\ l Ylm l (ux), specified by a global vec- 
tor ux=uiXi+U2X2- The ratio of u\ to u 2 characterizes 
the coordinate which is responsible for the rotation of the 
system [HJ HH. (The norm of u, uf+u^, simply affects 
the normalization of the basis function but not the rota- 
tion itself, and it may be set to unity. ) The isospin part of 
the system is expressed by t]tm t ■ The p-a Coulomb po- 
tential and the neutron-proton mass difference give rise 
to the isospin impurity in 6 Li, but its effect is rather 
small 1 1 81 ] . We ignore the neutron- proton mass differ- 
ence and consider no isospin mixing in the present study. 
The action of P\2 on the basis function is simple; it main- 
tains the functional form [l~5|, [lj| as follows: 

P 12 e-^ Ax [y L (ux) XS (l, 2)] JM vtm t (1, 2) 

= (-lf+ T e-^ A * [y L (tix) X3 (.l,2)}j M r,TM T M, 

(5) 

where the symbol ~ indicates the interchange of the el- 
ements of the matrix, that is, A\\=A22, Ai2=A2i=Ai2, 
A 2 2=An, and ui=u 2 , U2=u\. 

A conventional choice for describing the rotational mo- 
tion is a partial wave expansion, [y^ (xi)y^ 2 (x2)]lm l ■ 
One introduces a set of important partial waves (l -|l 2 ) 
to obtain a converged solution 0, [2(| ■ See also Ref. [L7| | 
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for the importance of including high partial waves in the 
Faddeev calculation for 6 Li. Our angular function looks 
quite different from this and it is much simpler than the 
partial wave expansion. No apparent couplings appear, 
and the use of the global vector greatly simplifies the 
calculation of the Hamiltonian matrix elements [H, • 
The merit of this angular function lies not only in its 
simplicity but also in its performance. As was compared 
in Ref . [16] , the global vector representation gives results 
as accurately as the partial wave expansion. Using the 
decomposition 

yLM L (uiXi+U2X 2 ) 
L 



4tt(2£ + 1)! 



e=o 



(2£+l)\(2L-2£+l)\ 



xu\u% t [yi(xx)y L -i(x 2 )]LM L , 



(6) 



the partial waves (£1(2) contained in the global vector 
part are only (OL), (1 L — 1), . . . , (£0). A very impor- 
tant point of our basis function is that other neces- 
sary partial wave contributions are brought about by 
the exp(— A\iX\ ■ x 2 ) term. When the term is ex- 
panded as J2 n [(— A 12 ) n /n\](xi • x 2 ) n , each term in the 
series produces the partial wave components of type, 
(xix 2 )"- A [3Ma;i)3Ma; 2 )]oo with A=n,n-2,...,1 or 0. 
When this is combined with [ye(xi)yf J -i(x 2 )]r J M L com- 
ing from the global vector part, it is clear that the ba- 
sis function (|U) can practically contain important partial 
wave combinations. 

We note that the basis function of Eq. (T3]) has a def- 
inite parity (— 1) L . As the ground states of 6 He and 
6 Li have a positive parity, this basis function cannot be 
used for L=l. We need to extend the basis function 
to make it possible to include L=l and a positive par- 
ity. This is made possible by replacing yLM L (ux) by 
[yL{ux)yi(u' x)]lm l [HI- For the case of two nucleons 
with L=l, this replacement results in a new basis func- 
tion 

V JM (A=(lS),A,x) = (I-P12) 

x {e-^iiy^xMx^xsMlMVTMAl^)}- 

(7) 

Note that [yi(xx)yi(x 2 )]i is equal to the vector product 
Xixx 2 within a constant factor. The operation of P\ 2 can 
be done as in Eq. ([5]), with an extra minus sign coming 
from this vector coupling. 

The set of A = (LS) included in the present calculation 
is summarized as follows: 



For 6 He (J 7 "" = 0+), 
For 6 Li (,r = 1+), 



(LS) = (00), (11). 
(LS) = (01), (10), (11) 



(21). 



(8) 



Here the basis function is given by Eq. (TjQ) for even L 
and by Eq. ([7|) for odd L, respectively. Note that the 
two independent L=l basis states exist for 6 Li. 



As noted before, the redundant state of the KKNN po- 
tential has to be eliminated; this elimination is a manifes- 
tation of the Pauli principle for the motion of the valence 
nucleons. The requirement is met by imposing that the 
trial wave function has no overlap with the 0s]y 2 bound 
state of the KKNN potential: 



(Qs; 



M*Jm)=0 (» = 1,2), 



(9) 



where the radial coordinate of the 0s]y 2 function is Xi and 
771 3 is the third component of the isospin of the nucleon 
i. The 0si/ 2 bound-state wave functions for the neutron 
and the proton are only slightly different because of the 
Coulomb potential, and thus we ignore their difference 
in the present calculation. The exclusion of the Osi/2 
component is practically achieved by the orthogonal pro- 
jection method .22!]. That is, we add in the three-body 
Hamiltonian |T]) a nonlocal, pseudo potential, 



i=i m=±i 



IOsj 



(i))(0s 



,(01. 



(10) 



and obtain a stable solution for sufficiently large 7 (typ- 
ically 10 6 MeV or larger than that). 

To search for good basis functions, we use the 
algorithm called the stochastic variational method 
(SVM) 15]. The SVM increases the basis dimension one 
by one by testing a number of candidates that are cho- 
sen randomly. Since each basis function is specified by 
the parameters (An, An, A 22 , m) or (An,Ai 2 ,A 22 ), the 
candidates are actually generated by giving random num- 
bers to the parameters chosen from physically important 
multi-dimensional parameter space. In this way we deter- 
mine about 100-200 basis functions for each A. The SVM 
works efficiently to take care of both the short-ranged re- 
pulsion of the realistic force and the elimination of the 
redundant states. 



C. Transformation of coordinate sets 

In the Faddeev method for a three-body system, three 
sets of T-type Jacobi coordinates are used, and each Fad- 
deev component of the total wave function, expressed in 
one of the three sets, is expanded in partial waves. To 
specify the basis function, however, we use the V-type co- 
ordinate x, which is different from the T-type coordinate. 
The V-type coordinate is a set of "single particle" (s. p. )- 
like coordinate, and it is chosen to make it easy to im- 
plement the symmetry of the nucleons. See Eq. (T5|). The 
s.p. coordinate is useful to represent the individual mo- 
tion of the nucleons when the correlation term vanishes, 
i.e., A\2 is set to zero. For example, the + ground state 
wave function of 6 He, approximated in p-shell harmonic- 
oscillator configurations, will be expressed in terms of a 
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combination of the spin-singlet and spin-triplet states 



D. Charge and matter radii 



$o =W e"5 oo(Xl+x5) 

x [[yi(x 1 )yi(a: 2 )]oXo(l,2)] 00 » 7 u(l > 2) J 

x [CM^iMiXiMlooWM), (11) 

where Ao and Afi are respective normalization constants. 
Here the shell model is extended to allow for different size 
parameters, ao and oi, for both the components. 

The T-type coordinate p—(r,R) is also convenient to 
impose the exchange symmetry for the two nucleons as 
Pi 2 simply changes r — > — r in the orbital function. (For a 
system consisting of more than three particles, the sym- 
metry requirement will be performed more easily in the 
V-type coordinate than other Jacobi coordinate sets.) 
Another set of coordinates commonly used is C=(Ci> C2)) 
called Y-type, which is related to the V-type coordinate 
by Ci=Xi and C2 =a; 2~ ^ where A c is the mass 

number of the core nucleus (A C =A in the present case). 
Each type of coordinate sets emphasizes its characteristic 
motion. As mentioned above, the V-type is suited for de- 
scribing the s.p. like motion around the core nucleus [23j ] . 
while the T-type coordinate is suitable for describing the 
motion corresponding to a+(2N) cluster decomposition. 
The Y-type coordinate plays a role similar to the V-type. 
In the limit of large A c , both the Y- and V-type coordi- 
nates coincide. 

It should be noted that the basis functions, ^ and 0, 
maintain their functional form under the transformation 
of the coordinates. The transformation from x to £, e.g., 
is done by a 2x2 matrix T as x=TQ. Then, it is easy to 
see that the basis functions change as follows: 

® JM (A,A,u,x) = 9 JM (A,f AT,fuX), 

*jm(A= (IS), A, x) = detr* JM (A= (IS), TAT, fl. 

(12) 

This flexibility of the basis function enables one to take 
account of possible important correlations of the system, 
just by choosing the nonlinear parameters suitably in 
only one particular coordinate set. 

Owing to the transformation property of Eq. (TT3]), the 
evaluation of the matrix element of an operator may be 
made in any convenient set of the coordinates. However, 
there is one exception. The matrix element for an angu- 
lar momentum-dependent operator such as the spin-orbit 
potential of Ui and the pseudo potential (flC)]) . if calcu- 
lated in the V-type coordinate, contains some error [24| . 
and it should be evaluated in the Y-type Jacobi coordi- 
nate set £. If the mass number of the core nucleus A c is 
sufficiently large as in our previous cases 25j, £2 is very 
well approximated by X2, and the error becomes small. 
In the present case, however, A c is only four, and we will 
see in Sec. IIII Al that the error is significant. 



The charge radii of 6 He and 6 Li can be calculated by 
taking into account the charge radii of the constituent 
particles of a and N as follows: 



,.2/6 



1 



( b Ke) = -(R 2 )+r 2 c (a)+r 2 c (n), 



9 

^Li)^ ) + 1 L 



r 2 ) + \r 2 c {a) + ^r 2 c {n) 



(13) 

\rlip)- 
(14) 




constituent particles are 
V^(p)=0.875fm, and 
Equation (fT5)) enables one 



T he ch arge radii 
v /rf(a)=1.671fm 
r 2 (ri)=-0.1161fm 2 
to deduce model-independent information on the (R 2 ) 
value from the experimental charge radius of 6 He. This 
will be discussed later. We can extend Eq. (fT5|) to a 
core+n+n system: 



(core+n+n) 



(R 2 ) + r 2 (core) 



Z 



!(n), 



(15) 



where A is the mass number of the system. In a 
9 Li+n+n three-body model for 11 Li, the recent data on 
the charge radii of 9,11 Li [28|, [2f| allow us to deduce 
VCR2y=5.95(3)fm. 

It is also interesting to calculate the point matter ra- 
dius which is defined by an rms radius for the point nu- 
cleon distribution. The point matter radii of 6 He and 6 Li 
are calculated by 

r^( 6 He, 6 Li) = 2 -(R 2 ) + l(r 2 ) + 2 -r 2 rn (a), (16) 

where r 2 n (a) is the mean square matter radius of the 
a particle and its value can be given by r^ i (a)=r~(a) — 
r 2 (p)—r 2 (n) under the isospin symmetry that the protons 
and the neutrons in the a particle have the same mean 
square radius. In what follows we use yj (a) = 1 .464 fm. 
Using Eqs. ([14")) and (fT!))) enables us to relate the matter 
radius of 6 Li to its charge radius: 



rf n ( b U)^r 2 c CU)~r 2 c (n)-ri(p). 



(17) 



The matter radius of 6 Li is thus expressed by only mea- 
surable quantities. 



E. Electric quadrupole moment 

The electric quadrupole moment operator is 
A 



16tt 



(18) 



where is the nucleon charge and X is the center of 
mass coordinate of the system. This operator can be 
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simplified, in the a+n+p model of 6 Li, to 



Q = Qr + Q R = \j^e\^y 20 {r) + ^ 20 (R) | . (19) 

Here use is made of the fact that the a particle has spin 
zero and the two valence nucleons have good isospin. The 
operator Q r is the same as the quadrupole moment op- 
erator for the deuteron, while the second term Qr corre- 
sponds to the quadrupole moment operator for the rela- 
tive motion between the a particle and the center of mass 
of the np system. 

F. Two-nucleon correlation function 



The density distribution for the TV- TV relative motion is 
given by the diagonal clement of the density matrix 

P(r) = g(r,r) = J W(r,k)dk, (23) 

and the momentum distribution for the N-N relative mo- 
tion is obtained by 

p(k) = J W(r,k)dr. (24) 

These distributions are normalized as J p(r)dr = 1, and 
J p(k)dk — 1. A formula to calculate the density matrix 
and the momentum distribution is given in Appendix. 



We define the two-nucleon correlation function by 

p(x 1 ,X 2 ) = 77J—7 ^{^Jm\xiX2){x 1 X 2 \^Jm)sT- 

M 

(20) 

Here (. . -) ST indicates that the integration is to be per- 
formed over the spin and isospin coordinates. Because of 
the average procedure over the z component of the to- 
tal angular momentum, the two-nucleon correlation func- 
tion becomes scalar, that is, it is a function of x\, x 2 and 
9, the angle between x% and x 2 ; p{x 1 ,x 2 )=p{xi,x 2l 9). 
The normalization of the total wave function ^ jm is ex- 
pressed as JJ p{x\,x 2 ,ff)%^ 2 x\x\ sin8dxidx 2 d9=l. 

G. Momentum distribution 

The momentum and spatial density distributions of 
a quantum-mechanical system are obtained through 
the Wigner phase-space distribution function [8f. The 
Wigner function is concisely expressed in terms of the 
density matrix. Since we are interested in the momentum 
distribution for the relative motion between the valence 
nucleons, we introduce the density matrix with respect 
to the relative distance vector r, one of the T-type coor- 
dinates: 

Q(r,r') = -±—J2 I (*JM\r'R)(rR\y JM ) ST dR, 

M 

(21) 

where (rR\ty jm) is obtained from (xix 2 ys> jm) through 

the replacement (fT2|) with T — ( \ ^ . We define the 

density matrix by taking the average over the z compo- 
nent of the total angular momentum. 

The Wigner distribution function is defined through 
the density matrix as 

W r(r 1 fc) = ^/ ff (r + f,r-f)e*-d.. (22) 



III. RESULTS AND DISCUSSIONS 

The input parameters we use are h 2 /mjv=41.47 
MeVfm 2 , m a =4mjv, and e 2 = 1.440 MeVfm, where mjv 
and m a are the masses of the nucleon and the a particle, 
respectively. The u parameter of the MN potential is set 
to u=l otherwise mentioned. No spin-orbit component 
of the MN potential is included. 

A. Spectroscopic properties 

In a single channel calculation with a specific A=(LS), 
neither (00) nor (11) channel produces a bound state of 
6 He below the a+n+n threshold. The case of 6 Li is dif- 
ferent, depending on the N-N potential. When using the 
MN potential, the A=(01) channel gives a bound state, 
whereas in the case of the G3RS potential no single- 
channel calculation produces a bound state below the 
three-body threshold. This difference between the two 
potential models is due to the tensor force which plays no 
role in the single (01) channel but gains energy through 
the coupling with different channels. This mechanism 
is similar to binding the neutron and the proton in the 
deuteron with the realistic force. 

Full calculations which couple all possible A channels 
of Eq. (JSj) give the results listed in Table HI The calcula- 
tion has been performed in the Y-type coordinates. The 
basis dimension K is 400 for 6 He and for 6 Li with the 
MN potential, while it is increased to 600 for 6 Li with 
the G3RS potential. The ground states of 6 He and 6 Li 
obtained with the G3RS potential are both underbound 
by 400-500 keV, while the MN potential underbinds 6 He 
by more than 500 keV but overbinds 6 Li by 200 keV. The 
underbinding of 6 He with the MN potential is partly due 
to the fact that the 1 So potential is too repulsive to re- 
produce the experimental phase shifts. The common 
lack of the binding energies in the case of the G3RS po- 
tential can be explained by at least three effects: One 
is the deficiency of the attraction in the D and F waves 
of the KKNN potential as discussed in Refs. 0, II]- 
According to the latter the energy gain of 6 He which is 
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TABLE I: Properties of the ground states of 6 He and 6 Li. Energy and length are given in MeV and fm. The L 2 and quadratic 
LS terms of the G3RS potential are neglected. See text for the MMN potential. The Coulomb potential is included in the 
term ( + U% ). Experimental data: £=-0.975, ^=2.054(14) [U for 6 He; £=-3.70, ^=2.55(4) [H[ or 2.540(30) [HJ, 
7^=2.42(4) for 6 Li. 







6 He (0+) 


6 Li 


(1 + ) 


d 


1+) 




MN 


MMN 


G3RS 


MN 


G3RS 


MN 


G3RS 


E 


-0.421 


-0.975 


-0.460 


-3.91 


-3.31 


-2.20 


-2.28 


(Tr) 


10.87 


11.87 


12.51 


17.56 


23.28 


10.48 


16.48 


<«8l> 


-3.77 


-4.86 


-5.62 


-13.41 


-7.71 


-12.69 


-7.29 




- 


_ 


0.107 


- 


-12.25 


- 


-11.46 




- 


_ 


0.021 


- 


- 


- 


- 


(Tr) 


12.47 


13.06 


12.55 


13.29 


11.49 


- 


- 


( u? + U§ ) 


-17.54 


-18.51 


-17.71 


-19.00 


-16.44 


- 


- 


{ Ut s + (7 2 LS ) 


-2.46 


-2.54 


-2.32 


-2.34 


-1.69 


- 


- 




5.05 


4.63 


4.86 


3.48 


3.58 


3.90 


3.96 




3.89 


3.66 


3.78 


3.49 


3.81 






V 'rri 


2.63 


2.49 


2.56 


2.27 


2.39 








2.09 


2.04 


2.07 


2.41 


2.52 






P(00) 


84.7 


86.4 


87.5 










p(ii) 


15.3 


13.6 


12.5 


1.1 


0.8 






P(10) 








6.2 


3.9 






P(01) 








91.7 


90.3 


100 


95.2 


P(21) 








1.0 


5.0 




4.8 



obtained by correcting the potential strength is, however, 
only a few tens of keV. Next is the effect of three-body 
forces [33[ though a conclusive statement on their mag- 
nitude remains open. The third effect to be considered 
is the distortion of the a core to 3iV + ./V partition or 
the clustering of the A=6 nuclei into 3 H+ 3 H (for 6 He) 
and 3 H+ 3 He (for 6 Li). The coupling of the a+N+N 
three-body configuration to the distorted configuration 
produces some energy gain. A recent microscopic calcu- 
lation indicates that the two configurations actually have 
rather large overlap and that the energy gain is of order 
of few hundreds keV in 6 He [HJ . 

We here remark on the accuracy of the present cal- 
culation by comparing to other calculations. For 6 He 
calculated with the MN potential, our ground state en- 
ergy, —0.421 MeV, is in excellent agreement with the 
value —0.42 MeV obtained by the hyperspherical coor- 
dinate method 35]. If the calculation is done in the V- 
type coordinates, (which is not correct as mentioned in 
Sec. Ill C|) . the ground state energy of 6 He would go down 
to —0.749 MeV; the result is again consistent with the 
values obtained by a Lagrange-mesh calculation [24| and 
a hybrid T+V model calculation 14j. The correspond- 



ing energy for the 6 Li ground state would be —4.68 MeV, 
instead of —3.91 MeV calculated correctly in the Y-type 
coordinates. These results clearly show that the correct 
treatment of the coordinates is important in such light 
systems as 6 He and 6 Li. 



Table |T] lists the decomposition of the energy into the 
expectation values of the various pieces of the Hamilto- 
nian, the radii and the probability P(LS) in % of finding 
the A=(LS) component. Let E r and Er denote respec- 
tively the energy for the relative motion of the two nucle- 
ons and the energy for the relative motion between the 
a particle and the center of mass of the two nucleons. 
They are defined by 



E r — \T r ) 
E R ={T R 



"12 / 



U 



LS \ 
; 12 / I 



u. 



LS 



LS 



(25) 



In the case of 6 He, the G3RS and MN potentials give 
very similar results for these energies: (E r , Er)=(7.03, 
-7.49) MeV for G3RS and (7.10, -7.52) MeV for MN. 
Even each of the expectation values is close to each other 
as well. The tensor and spin-orbit forces between the two 
neutrons play a minor role. In contrast to 6 He, both the 
potentials exhibit quite different features in binding 6 Li. 
The partial energies are (E r , E R )=(3.32, -6.63) MeV for 
G3RS and (4.14, -8.05) MeV for MN, so that the differ- 
ence between the two is modest. However, the content 
of E r in particular is quite different between them be- 
cause of the tensor force and the short-ranged repulsion. 
In G3RS the large positive value of ( T r ) is balanced by 



the tensor contribution 



Though this is similar to 



the case of the deuteron, we see from the values of ( T r ) 
and \J (r 2 ) that the np pair in 6 Li is more compressed 
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than the one in the deuteron. In spite of these differ- 
ences the P(LS) values of 6 Li are rather similar in the 
two potentials except for the (LS)= (21) channel which 
is largely determined through the tensor coupling to the 
dominating channel of (LS)— (01). As will be discussed 
later, the 6 Li quadrupole moment is very sensitive to this 
coupling. The spin-triplet channels occupy about 95%. 

The 6 Li charge radius calculated using the G3RS po- 
tential agrees with the experimental value determined 
from the electron scattering [3l[ and the isotope shift in 
lithium but the MN potential gives the charge ra- 
dius which is small by at least 0.1 fm. T his failure of the 
MN potential is due to that the \J (R 2 ) value is predicted 
to be too small. As Eq. (fT?)) shows, the use of the 6 Li 
charge radius leads to an "experimental" rms matter ra- 
dius \J r 2 n ( 6 Li)=2.42(4) fm, which is consistent with that 
determined from the analysi s of prot on elastic scatterings 
at intermediate energies, y / r£ l ( 6 Li)=2.45(7) fm [361 l37j|. 
The theoretical value calculated with the G3RS potential 
is 2.39 fm, in agreement with experiment. 

Now we turn to the case of 6 He. The charge ra- 
dius of 6 He has recently been measured by laser spec- 
troscopy technique [30]. Our theoretical value agrees 
fairly well with the experimental value. Using the ob- 
served charge radius of 6 He in Eq. I|13p. we can deduce 
the "experimental" value of yj (R 2 ) , which turns out to 
be 3.726(69) fm. Our value of 3.78 fm is very consistent 
with this value considering the fact that the calculated 
binding energy is small by 500 keV. The rms matter ra- 
dius of^Heextracted from the proton elastic scattering 
is \/r^( 6 He)=2.30(7) fm 37|, which is, however, smaller 



than the value (2.48±0.03fm) deduced from the inter- 
action cross section analysis [38j . The theory with the 
G3RS or MN potential predicts too large matter radius 
for 6 He; this overestimation is related to that the calcu- 
lated bindi ng en ergy is too small, leading to a too large 
value for y/ (r 2 ). Since the attraction of the MN potential 
is weak in the 1 Sq channel, we repeated the calculation 
by increasing the strength of its longest range part from 
-91.85 to -91.85x1.07 MeV so as to reproduce the bind- 
ing energy. This potential is called MMN, and its result 
is listed in Table Q] As expected, the rms matter radius 
now decreases from 2.63 to 2.49 fm, which is consistent 
with the values of the interaction cross section measure- 
ment as well as the fully microscopic three-cluster calcu- 
lations [ii m. 

The electric quadrupole moment of 6 Li is a long- 
standing problem. The 6 Li quadrupole moment is neg- 
ative and small, -0.08178(164) efm 2 [39[, while the 
deuteron quadrupole moment is 0.2860(15) efm 2 [io| . 
Thus it is known that the a+n+p model which assumes 
(0s) 4 a-cluster gives the wrong sign even when the tensor 
force is included @. See also Ref. [4l|. Though we do 
not attempt to solve this enigmatic issue, we just show 
the result of the present model in Table [TTJ As expected, 
the 6 Li quadrupole moment obtained with the G3RS po- 
tential turns out to be positive (0.164 efm 2 ). Interest- 
ingly, the quadrupole moment with the MN potential 



TABLE II: Electric quadrupole moments of 6 Li and the 
deuteron in efm 2 . Values in the parentheses are the con- 
tributions of the cross terms between the (LS) = (01) and 
(LS) — (21) components of the 6 Li ground state wave func- 
tion. The experimental values are —0.08178(164) efm 2 [3^] 
for 6 Li and 0.2860(15) e fm 2 [U for the deuteron. 





6 Li (1+) 


d(l+) 




MN G3RS 


G3RS 


Q 


-0.295 0.164 
(-0.396) (0.088) 


0.264 




-0.034 0.198 
(-0.094) (0.160) 


0.264 


Qr 


-0.260 -0.034 
(-0.302) (-0.072) 





becomes —0.295 efm 2 ; this happens because the large 
negative contribution of Qr is not canceled by the Q r 
value. This contrast between the two potential models is 
mainly due to the different contribution of the cross terms 
between the (LS) = (01) and (LS) — (21) components of 
the 6 Li ground state wave function. See the parenthetic 
values in Table II. The mixing of these components is 
due to the tensor force. Therefore the failure of the 6 Li 
quadrupole moment indicates that we have to consider 
the effect of the tensor force between the core and the 
valence nucleons. It should be noted that a variational 
Monte Carlo calculation gives the quadrupole moment of 
-0.23(9) efm 2 



B. Two-nucleon correlation function 

Figure [1] plots the density distributions p(r) (normal- 
ized to unity) of the two-nucleon relative motion in 6 Hc, 
6 Li and the deuteron. The densities calculated using the 
G3RS potential (right panel) show central dips due to 
the short-ranged repulsion, but beyond r=1.5fm they 
are similar to those calculated with the MN potential 
(left panel). The density of 6 He reaches furthest in the 
distance, and as a result its density around r=l~2fxn is 
considerably smaller than that of 6 Li. Comparing the 
densities between 6 Li and the deuteron, we see that the 
np relative motion in 6 Li shrinks compared to that of the 
deuteron (see also the y/ (r 2 ) value in Table EJ). 

The correlated motion of the valence nucleons reflects 
on the two-nucleon correlation function p(xi,X2,0). Be- 
fore discussing the correlated motion, we first examine 
the function p(x\, X2, 0) which is generated from an "un- 
corrected" basis function <!> for 6 He. For this purpose we 
take a combination of the two p-shell harmonic-oscillator 
functions (fTT|) . 



$ = \Zl-C 2 <f> + C<f> u (26) 
and determine (ao, ai, C) so as to maximize the overlap, 
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| ($|Woo) | 2 , with the 6 He ground-state wave function <3>oo 
obtained using the G3RS potential. The resulting values 
are a =0.163fm- 2 , ai=0.194f m - 2 , and C=0.402, lead- 
ing to | (<E>|\I/oo)| 2 =0.75. The simple wave function $ has a 
surprisingly large overlap with the realistic wave function 
^oo- Though the overlap is fairly large, $ includes no cor- 
related configurations and indeed the energy calculated 
with $ is high (8.77 MeV). The two-nucleon correlation 
function p(x\ 1 X2,Q) calculated from $ becomes a func- 
tion of cos 2 9, so that p(x, x, 6) multiplied by 87r 2 a; 4 sin# 
is symmetric with respect to 0=90°. See Fig. [2] An asym- 
metry with respect to 0=90° would indicate the presence 
of correlation in the A=6 nuclei. 

Now let us discuss the two-nucleon correlation func- 
tion which derives from the dynamical calculation with 
the G3RS potential. The MN potential gives similar two- 
nucleon correlation function. Figure [3] displays the con- 
tour maps of 8ir 2 x 4 smO p(x, x, 8) for 6 He and 6 Li. In 
both the cases we clearly see asymmetric patterns with 
two distinct peaks: In 6 He the highest peak is located 
at (£,0)=(2.9,26°) with a height of 0.07 fin" 2 , while in 
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FIG. 2: Contour map of the two-nucleon correlation func- 
tion, 8n 2 x 4 sin 8p(x, x, 8), calculated from the uncorrelated p- 
shell wave function, Eq. (|26[) . The difference between any two 
neighboring contour levels is 0.01 fm~ 2 . 



6 Li the highest peak is located at (x, 0)=(3.2, 24°) with 
a height of 0.13fm~ 2 . Here x is given in fm. The peak 
in 6 Li called the deuteron-like correlation is about twice 
higher than that called the dineutron-like correlation in 
6 He. Comparing Fig. [3] with Fig. [2j we learn that the 
two-nucleon interaction enhances the asymmetric pat- 
tern. Another lower peak which shows up at larger angles 
is called a cigar-like correlation. It corresponds to the ge- 
ometry that the two nucleons sit on the opposite sides of 
the core. This type of peak is located at (x, 6>)=(2.3, 
136°) for 6 He and at (x,6)=(2A, 136°) for 6 Li, respec- 
tively. Their heights are both 0.03 fm~ 2 , which is about 
half of the peak height in Fig. [2] 



C. Projection to dineutron- and cigar-like 
configurations 

In Sec. IIIIBI we discussed the two-nucleon correlated 
motion from the viewpoint of the asymmetric appear- 
ance of the dineutron- and cigar-like peaks. Here we ask 
a question of how the two peaks appear in relation to 
the density or rms distance distribution of the two nu- 
cleons. We expect that the dineutron-like peak in 6 He 
is formed from those components of the wave function 
which have smaller rms distance, whereas the cigar-like 
peak is constructed from the rest of the components. 

As Eq. ([3]) shows, the wave function is given as a combi- 
nation of K independent, nonorthogonal basis functions. 
By taking a suitable linear transformation of these K 
bases, we can obtain an orthonormal set Sjm(m), 

k 

2jm(m) = Y, W ^JM(KA,Ui), (27) 

i=l 

with the condition (SjM(A t )I^JM(A < ')) = ^^'- ^ ^ s con " 
venient for the present purpose to choose the coefficients 
in such a way that they diagonalize the squared dis- 
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FIG. 3: Two-nucleon correlation functions, 8tt 2 x 4 sin 9p(x, x, 
panels are their contour maps, and the difference between any two neighboring contour levels is 0.01 fnT 



calculated with the G3RS potential for 6 He and 6 Li. The lower 



tance r 2 between the two nucleons, that is, 
(Ej M {t i )\r 3 \Sj M (j i t )) = (r\5 llj , 



(28) 



r'l m increasing order 



By arranging the eigenvalues l^ 2 ^ 

(/i=l,2,. . . , K), we display in Fig. [4] the distribution of 
the rms distance in the case of 6 He, where K=AQQ and 
the minimum and maximum eigenvalues of W (r 2 ) M are 
0.223 and 40.7 fm, respectively. The fact that the eigen- 
values cover the wide region from small to large distances 
indicates that the SVM basis selection is efficiently per- 
formed to take into account the short-ranged correlation 
as well as the asymptotic behavior. We can see that the 
basis states S jm (h) of the first 320 members give rather 
uniform distribution of yf (r 2 ) p up to about 12 fm, while 
the rest of the basis states cover the eigenvalues of larger 
rms distances. 

We define two projectors which are orthogonal com- 
plements to each other: 



P S = J2 ISjm(m)) (E jm (h)\, 

K 

Pl= |Hj M (/i)) <3 JM (/i)| , (29) 



spanned by the basis states with larger rms values. The 
total wave function is decomposed into two orthogonal 
components, "small" and "large" , using the projectors 



JM 



JM 



Pj^ 



JM 



<F S + tf L . (30) 



For the sake of simplicity, we choose k in Eq. 
such that K^sI^jm)! is as close as to 0.5. It turns 
out that k=88 and the corresponding rms eigenvalue is 
v / (r2)^=3.47fm. 

It is now possible to decompose the expectation value 
of an operator O into three terms, that is, small, large 
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Label of basis \l 

FIG. 4: Eigenvalues of the rms distance between the two neu- 



with Ps+Pl=1. The projector Pg projects into the sub- ^5°^ s ' 
space spanned by those Sjm(m) which have smaller (r 2 )^ 
values, while Pl is the projector into the rest which is 



■ 2 ), calculated from the basis functions for He. The 
G3RS potential is used. Dotted line denotes the rms distance 
(4.86 fm) of 6 He. 
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and their interference terms: 

{^jm\0\*jm) = {*s\0\9 s ) + 

+ {(*s\0\*l) + (*l\0\V s )}. (31) 

We apply this decomposition to the two-nucleon correla- 
tion function to see how the contour map of 6 He (Fig. [3]) 
is constructed. Plotted in Fig. [5] are those contributions 
to the contour map which are calculated with tyg and \& l, 
respectively. The contribution of the interference term is 
found to be small and can be safely ignored. Compar- 
ing this decomposition with the full contour map of 6 He 
in Fig. [21 we can safely conclude that the dineutron-likc 
correlation is generated by the small component \& g and 
the cigar-like correlation by the large component it s- 



D. Momentum distribution 

Though the contour map discussed in Sees. IIII Bl and 
MI CI shows some correlation effects, it is not clear how 
the correlated features in 6 He and 6 Li are observed exper- 
imentally. Comparative experiments of the intermediate 
energy proton elastic scatterings on 6 He and 6 Li (3(| |37| 
have been performed in order to elucidate the matter den- 
sities of both the nuclei, but the analysis of the experi- 
mental data is confronted with some ambiguities because 
the scattering is confined to extremely forward angles. As 
mentioned in the Introduction, the measurement of the 
momentum distribution in a special arrangement seems 
to be accessible in the inverse kinematics, providing data 
which are sensitive to the different structures of 6 He and 
6 Li. 

It is well-known that the momentum distribution of the 
np relative motion in the deuteron shows different behav- 
ior in the S- and D-wave contributions. As displayed in 
the right panel of Fig. [6l the S'-wave contribution to the 
momentum distribution is peaked at lower momentum 
and has a dip at fc~2fm . The D-w&ve component of 
the deuteron, however, fills the dip. This characteristics 
of the distribution is supported by experiment (43[. In 
contrast to this, the momentum distribution (left panel) 
obtained with the MN potential does show a dip because 
it has no D-w&ve component, and in addition the mo- 
mentum distribution decreases quickly with increasing k 
because the short-ranged repulsion is not as strong as 
the G3RS potential. To compare with experiment at k 
higher than 2fm _1 , however, it is important to include 
meson exchange currents and isobar currents dominated 
by the A excitation. See Ref. and Fig. 

The momentum distributions of 6 He, 6 Li and the 
deuteron are compared in Fig. [7] for the G3RS (right 
panel) and MN (left panel) potentials. The realistic po- 
tential of G3RS gives the momentum distributions char- 
acterized as follows: The momentum distribution of 6 Li 
is very similar to that of the deuteron, but the momen- 
tum distribution of 6 He differs from them, showing a 
clear dip at A;~2fm _1 . These features are understood 



from the difference in the partial wave contents of the 
N-N relative motion; 6 Li contains the D-wave compo- 
nent as the deuteron does, whereas 6 He is dominated by 
the S'-wave component. The most distinctive difference 
between 6 He and 6 Li appears around fc~2fm . In this 
region, however, the momentum distribution becomes by 
four or more order of magnitude smaller than that at 
0; this may make it hard to measure the cross section 
experimentally. If the measurement of the momentum 
distribution is possible in this region, one can learn the 
role of the tensor force acting between the valence nucle- 
ons, provided that the meson exchange currents and the 
isobar excitations are still not so important. 

Figure [5] compares the momentum distributions of 6 Hc 
corresponding to the three different wave functions, those 
obtained with G3RS and MN, and the uncorrelated one 
defined in Eq. J5B]). Both the G3RS and MN distribu- 
tions are similar up to the dip region. Beyond fc^2fm _1 
the momentum distribution of G3RS surpasses that of 
MN, which is due to the difference in the short-range 
correlation involved in the two wave functions. The un- 
correlated wave function gives the momentum distribu- 
tion which is quite different from those of the correlated 
wave functions even at fc^ lfm -1 . 

In Sec. IIII CI we decomposed the ground state wave 
function of 6 He into ^>g and ^l, and confirmed that the 
dineutron-like peak is produced by the small component 
^5, while the cigar-like peak by the large component ^>l- 
The interference term was small. One might expect that 
the momentum distribution may as well be decomposed 
into low and high momentum components. Namely, 
may contribute to the momentum distribution at small k, 
while ^ g to the high momentum component. We exam- 
ine this expectation in Fig.[9]by analyzing the momentum 
distribution of 6 He obtained with the G3RS potential. 
The left panel shows the partial momentum distributions, 
PL calculated with and ps calculated with ^s, re- 
spectively. The right panel compares the full momentum 
distribution p with the incoherent sum of pl+Ps, so that 
the difference between p and pl+Ps is the contribution 
of the interference terms of Eq. (|3"Tj) . We see that the 
contribution of the interference terms can be neglected 
for k <lfm . However, the interference contribution 
becomes important for fc>lfm _1 . In the momentum re- 
gion where the interference can be ignored, the momen- 
tum distribution is dominated by for k < 0.5 fm -1 
and by $g for 0.5 < k < 1 fm -1 . 



IV. CONCLUSIONS 

To study the correlation and the momentum distribu- 
tion of the two-nucleon relative motion in the ground 
states of 6 He and 6 Li, we have described these states in 
a three-body model of a+N+N where the a particle is 
assumed to be an inert core. We used a parity-dependent 
a-N potential which reproduces the low-energy S- and 
P-wave phase shifts, and two different types of N-N in- 
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FIG. 5: Decomposition of the two-nucleon correlation function of 6 He into the small (left panel) and large (right panel) 
components. The G3RS potential is used. The difference between any two neighboring contour levels is 0.01 fm~ 2 . 
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FIG. 6: Momentum distributions of the deuteron calculated with the G3RS and MN potentials. Data are taken from Ref. [44 



teractions as the potential acting between the two va- 
lence nucleons. One is a realistic potential which con- 
tains the tensor and spin-orbit forces and the other is an 
effective potential which includes no tensor component. 
These were used to compare how much the different N- 
N potentials affect the correlation and the momentum 
distribution. 

We have obtained the solution of the three-body prob- 
lem by approximating the 6 He and 6 Li ground state wave 
functions in terms of a combination of explicitly corre- 
lated Gaussian basis functions. The use of the global 
vectors to describe a nonspherical orbital motion facili- 
tates the calculation of the matrix elements much easier 
than the partial wave expansion, and moreover provides 
us with a solution of high accuracy. 

The energies and rms radii of 6 Hc and 6 Li are com- 
pared to experiment. The energies calculated with the 
realistic N-N potential are underbound by 400-500 keV 
in both the cases. The charge radii of 6 He and 6 Li and 
the matter radius of 6 Li are in fair agreement with the 
observed values. The 6 He matter radius is predicted to 



be larger than experiment; the result is due to the under- 
binding of the calculated ground state of 6 He. We have 
analyzed the two-neutron correlated motion in 6 He in 
order to identify how the dineutron and cigar-like config- 
urations are related to the two-neutron relative distance 
distribution. 

The momentum distributions of the N-N relative mo- 
tion have been compared between 6 He and 6 Li. The 
distributions obtained with the effective potential show 
the pattern characteristic of S-wave dominance and fall 
rapidly as the momentum increases. In the case of the 
realistic potential, the momentum distribution in 6 Li is 
very similar to that of the deuteron. That is, both the S- 
and D-waves contribute to the momentum distribution 
which monotonically decreases with an increasing mo- 
mentum. In contrast to this, the 6 He momentum distri- 
bution is dominated by the S*-wave, showing a clear dip at 
fc~2fm . The most prominent difference in their mo- 
mentum distributions thus shows up around fc=2fm _1 . 
The difference between 6 He and 6 Li is primarily due to 
whether or not the tensor force plays an important role 
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FIG. 7: Momentum distributions, calculated from the G3RS and MN potentials, of the valence nucleons in 6 He and 6 Li. 




FIG. 8: Momentum distributions of the valence nucleons in 
6 He for the three different wave functions. 



of mixing the Z?-state probability between the N-N rela- 
tive motion. We hope that this prediction will be tested 
experiment ally. 
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APPENDIX: CALCULATION OF THE 
MOMENTUM DISTRIBUTION 

The aim of this appendix is to show a method of cal- 
culating the momentum distribution for the correlated 
Gaussians. As discussed in Sec. Ill Gl the momentum dis- 



tribution is calculated from the density matrix, so that 
it is sufficient to show how the density matrix is evalu- 
ated. We first express the basis functions in the T-type 
coordinates (zi,z 2 ) as discussed in Sec. Ill CI (in this ap- 
pendix we use (21,2:2) instead of (r,R) to simplify the 
notation), and write the general form of the orbital part 
of the correlated Gaussians as 



G L m l {A, £1, £2, zi, Z2) 
= exp (-^zAzj [yi 1 {zi)y l2 {z 2 )] LML 



(32) 



The density matrix p which we consider here reads as 
P(zi,z[) 

= E ^ Gl ' e » *a)xs'(l,2)] JM | 

/J + 1 M 

x \[G L (A,£ 1 J 2 ,z 1 ,z 2 )xs(l,2)]j M ), (33) 



where the integration in the matrix element has to be 
performed for 22 as well as the spin coordinates. 

Performing the integration over the spin coordinates 
yields 



p{z\,z[) = 5 LL 'Sss' 2L , 1 E 



x (G LMl (A', £[, z{, z 2 )\G LMl (A, l x ,£i, z 1 ,z 2 )) 



(34) 



Writing the angular part of the right-hand side of Eq. 
explicitly we obtain 

^— £ [yMy i2 { Z2 )] LML [y q (z{)y q (z 2 )Y LML 

Ml 

= ( _l)4+4_^ = y 
V2i + 1 ^ 



V2L 

[y,Azi)yq(z{)] x yx(z2 



£1 £2 L 
£{ V 2 L 
A A 



C(£ 2 £^X)4 



on 



(35) 
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FIG. 9: Momentum distribution of the valence nucleons in 6 He and its decomposition into small and large components. See 
the text for detail. The G3RS potential is used. 



where [ ] is the 9j symbol in unitary form and C is the 
coefficient to couple two spherical harmonics with a same 
argument: 



Note that 2n = £ 2 +£ 2 — ^ i s non-negative and even, oth- 
erwise the coefficient C(£ 2 £ 2 ; ^) vanishes. Thus the inte- 
gration over z 2 in Eq. (|34[) is performed as 



exp 



1 



-az 



Z-z 2 z 



[y tl {zi)y l{ {z{)] x y x {z 2 ) 



dz?_ 



(10 



tt (2n)!! (A+i) 



00 



1 exp , 
2a J V 2a 



(37) 



where L„ 2 is the associated Laguerre polynomial, and 



a = A 22 +A!, 2 , Z = A 12 zi + A{ 2 z[. 



(38) 



Using the formula $Q the coupling of three y's in 
Eq. (j37|) is reduced to 



[yx(z) [ytMytM)]. 

= E 



110 



4tt(2A + 1) 



(2fc+l)!(2A-fc + l)! 



4 fc 4 ' 



12^12 



E 



A— A; A 
4 # A 



(i [i 



C(k£v,Li)C(\-k£l;iJ,) 



(39) 



with 



MzOi^OU = (-ir 



47T 



K=0 



(2/x-2/t-l)!! /zrZj' 



(/i-2«;)!(2«;)!! V ^ 



(40) 



where [#] is the largest integer less than or equal to ^. 

Combining Eqs. (fM| - (F4"0"|) , we obtain the density ma- 
trix as a combination of terms 



zl p zi 2p ' ( Zl -z{y exp (~/3z 2 - fi'z? - yz v zi) , (41) 



where p, p' and g are all non-negative integers. 

To calculate the momentum distribution we just re- 
place (zi,z{) with (r + i«,r — |s) in the density matrix, 
multiply e ifc s and integrate over r and s. Renaming 
(r, s) as (zi, z 2 ) again, the integration results in the fol- 
lowing form 



e-i SBz+ii * z zf ni z 2 2n2 {zi-Z2) n3 dzidz2, (42) 



which can be performed analytically, where B is a 2 x 2 
symmetric matrix, kz —h-Z\ — h-z 2 and n\, n 2 and 713 
are all non-negative integers. 
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